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The torque in turbulent Taylor- Couette flows for shear Reynolds numbers Res up to 
3 x 10 4 at various mean rotations is studied by means of direct numerical simulations 
for a radius ratio of r\ = 0.71. Convergence of simulations is tested using three criteria of 
which the agreement of dissipation values estimated from the torque and from the volume 
dissipation rate turns out to be most demanding. We evaluate the influence of Taylor 
vortex heights on the torque for a stationary outer cylinder and select a value of the aspect 
ratio of T = 2, close to the torque maximum. The connection between the torque and 
the transverse current J u of azimuthal motion which can be computed from the velocity 
field enables us to investigate the local transport resulting in the torque. The typical 
spatial distribution of the individual convective and viscous contributions to the local 
current is analysed for a turbulent flow case. To characterise the turbulent statistics of the 
transport, PDF's of local current fluctuations are compared to experimental wall shear 
stress measurements. PDF's of instantaneous torques reveal a fluctuation enhancement 
in the outer region for strong counter-rotation. Moreover, we find for simulations realising 
the same shear Re$ 2 x 10 4 the formation of a torque maximum for moderate counter- 
rotation with angular velocities lu sa — 0.4 Wj. In contrast, for Res ^ 4 x 10 3 the torque 
features a maximum for a stationary outer cylinder. In addition, the effective torque 
scaling exponent is shown to also depend on the mean rotation state. Finally, we evaluate 
a close connection between boundary-layer thicknesses and the torque. 
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1. Introduction 

The flow between two concentric ind ependently ro tating cylinders shows, after a first 
bifurcation to axially periodic vortices (|Tavlorlll923l ). various routes to turbulence. The 
presence of two easily accessible control parameters, the rotation rates of the two cylin- 
ders, and two geometrical parameters, the ratio of the radii and the aspect ratio, the 
height of the cyl inders relative to the gap, provide a , huge parameter space with a rich 
phenomenology ( Andereck et al. 1986t Koschmieder 1993 : rChossat &; Ioossl 1994 ). For 
the turbulent situation one might expect much less of a dependence on these parame- 
ters because of the homogenising effects of turbulent fluctuations, but already the ex- 
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periments bv IWendd (1 19331). iTavlorl (Il936l ) and later studies dLathrop et all Il992^f6| : 
Lewis fe Swinnevlll999l: iRavelet 6^120101 : iBurin et a^2010l : IPaoletti fc Lathropll201ll : 



van Gils et a/.ll2011g| |6l) show that there is considerable structure in the data that needs 



to be explained. 

The measurements of torque bv IWendU ( 19331 ) up to Reynolds numbers 10 5 show that 
for the case of the outer cylinder at rest the dependence of the torque on the inner cylinder 
velocity can b e described by two effective power laws, one for high and one for low shear 
rates, see also ( Racina fc Kindf2006t Burin et aLll2010l ). In contrast, high-precision torque 
measurements with a stationary outer cylinder revealed that the torque does not follow a 
pure power-law scaling since the local scaling ex ponent increases monot onically with the 
inner cylinder velocity (jLathrop et al. 1992 alia: iLewis fc Swinnevlll999l) . Moreover, first 
torque determinations also for counter-rotating cylinders and for the inner cylinder at 
rest showed a n even more complic ated dependence of the torque scaling exponent on the 
applied shear (jRavelet et «Lll201(J ). In addition to the scaling with the shear rate, recent 
experiments for Reynolds numbers up to a few 10 6 revealed that, when the mean shear 
rate is kept constant, the maximum in torque i s reached for counter-rotating cylinders 
(|Paoletti fc Lathrodl201ll Ivan Gils et aZj^Olld fih. 

These high Reynolds number experiments were complemented in narrow Reynolds 
number regimes with direct numerical simulations (DNS) of turbu lent Taylor-Couette: 
Flow characteristics were analysed for the outer cylinder at rest ( Bilson fc Bremhorsd 
20071 : lDongl2007t IPirro fc Quadriol2008l) . for counter-rotating cylinders (lDond20086l) an d 



in the case of spiral turbulence (M eseguer et alJl2009t iDon dl2009t iDong fc Zhendl201lh . 
The highest Reynolds number (~ 10 4 ) was achieved bv IPirro fc Quadriol (|2008j T and it 
touches the Reynolds number range studied experimentally, but they do not analyse the 
turbulent transport and its contributions to the torque. Our aim here is to present direct 
numerical simulations of Taylor-Couette flow up to Re = 3 x 10 4 and to show that with 
these simulations one can explore the cross-over regime from the low, vortex dominated 
Re regime to the more fully turbulen t regime accessible in experi men ts. We will be guided 
by previous numerical simulation s oflCoughlin fc Marcusl (Il996l) andlDona (120071 12 008b) 



and the experimental findings of Paoletti fc Lathropl ( 2011 ) and van Gils et al. ( 2011al ). 

The outline of the paper is as follows. In section [2] we define the problem, describe 
the numerical method and explain the convergence tests we applied. The dependence 
of the torque on the domain height is analysed in section [3] The local contributions 
to the torque near the inner and outer cylinder and the variations with the radius are 
discussed in s ection |H This sect i on als o includes a discussion of the fluctuation statistics 
measured bv ILewis fc Swinnev ( 19991 ). In section [5] we discuss the dependence of the 
torque on bot h Reyno l ds numbers and dis c uss th e relation to rec e nt ex periments by 
van Gils etaJ\ (|2011al ): IPaoletti fc Lathropl (|201lh : Irluisman et all (|2012h . Finally, we 



discuss the determination of a bou ndary-layer thickness th at is compatible with the 
definition of the current proposed in lEckhardt et all ( 20076 ) in section [6] We conclude 
with a brief summary and outlook. 



2. Problem setting and numerical methods 

2.1. Definitions 

We investigate the incompressible flow between two concentric cylinders by means of 
direct numerical simulations. The geometry of the Taylor-Couette (TC) system is char- 
acterised by the radius ratio n = ri/r Q , where r, and r denote the radius of the inner 
and outer cylinder, respectively. Here, we study the case T] — 0.71 which is close to values 
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chosen in recent experimental torque studies ([Lathrop et al. 1992 ol fi!: iLewis fc Swinnev 



19991 : Ivan Gils et aLll2011a : Paoletti fc Lathropll201ll ) The simulations are performed in 



an axially periodic domain and therefore do not include the effects of top and bottom lids 
present in TC experiments. Unless otherwise specified, the aspect ratio T — L z /d = 2 is 
chosen, with L z the height of the cylinders and d the gap width, d = r D — ri. This domain 
can support one Taylor vortex pair when the outer cylinder is at rest. In the figures we 
employ the rescaled radial coordinate y = (r — r^/d which ranges between for the inner 
cylinder and 1 for the outer cylinder. 

For the measure of the shea r between the cylinde rs we use the definition of a shear 



Reynolds number Res given in lDubrulle et al\ ([20051 ) . 



Re s = -—\r ] Re -Re i \, (2.1) 

1 + 7? 

with Rei = driUJi/v and Re Q = dr uj jv^ where oji and to denote the angular velocity of 
the inner and outer cylinder and v the kinematic viscosity of the fluid. We measure the 



mean (or global) rotation of the system by the ratio of angular velocities ([van Gils et 
201 lab 



a = —ujo/uji. (2-2) 

Note that for Re = the shear Reynolds number Res — 2/(1 + rf)Rei is slightly larger 
than Rei. 

To connect to the analysis in lEckhardt efal\ (|20076l) . we consider the local current of 



angular velocity 

j"{r,<p,z,t) =r 3 (u r u~isd r Lj), (2.3) 

where u r and u v denote the radial and azimuthal velocity components and w = u v /r is 
the angular velocity. Its average over any cylindrical surface concentric with the bounding 
cylinders gives the mean angular velocity current 

J" = U") A{ r),t = r 3 ({uM A{r)>t - vd r (w) A(r)it ) , (2.4) 

which is related to the dimensionless torque G — T / (2nL z pfV 2 ) — ^~ 2 J W , where T is 
the torque needed to drive the cylinders and pf is the density of the fluid. Here (• • -) A / r \ t 
stands for averages over a cylindrical surface of radius r within ^ r ^ r Q and over 
time. As the current is independent of r, any ^ r $J r Q can be used to calculate it, but 
the rate of convergence varies with r (see below) . 

A rescaled value of the angular velocity current which is analogous to the Nus selt 



number describing the heat current in thermal convection (|Eckhardt et al\ l2007af ) is 
obtained by dividing J u by its laminar value for the circular Couette flow, Jf am = 
2vr?r 2 (u)i-u )/{rl-r?), 

Nu u = .r/J? am = G/G lam . (2.5) 

The resulting quantity Nu^ is the quasi-Nusselt number and gives the torque in units of 
the laminar value. 

In the following we will study the contributions of the different terms in (|2.4|) to the 
local and global values, the fluctuations of the local currents and the variation of the 
global torque with Reynolds number and rotation ratios for shear Reynolds numbers up 
to 3 x 10 4 . 

2.2. Numerical considerations 
The n umerical code used for all subsequent computations was developed by lMeseguer et al. 



(|2007l ) and employs a solenoidal spectral basis to accomplish a Petrov-Galerkin scheme, 
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e.g. Moser et al. ( 19831 ): Meseguer ( 2003 ). For the computations all variables and fields 



are rendered dimensionless with the gap width d and the viscous time d 2 jv as character- 
istic scales for space and time. In the numerical scheme the velocity field u = m& + v is 
decomposed into its laminar contribution Ub and a deviation v. The latter is expanded 
in a Fourier- Chebyshev vector field basis $; nm , 



L N M 



v{r,ip,z,t) = V V Ya lnm (t)$ lnm , &i nm = e i V koZ+nv) v lnm (r), (2j 



l=-L n=—N m=0 



with ko = 2n/T and spectral coefficients a;,, 
satisfy the incompressibility condition V • <&z nm = and boundary condition vi nm (ri) = 
vinm{fo) — and the vector functions Vi nm are constructed from Chebyshev polynomials. 
The highest order of modes employed in axial, azimuthal, and radial direction is given 
by L, N, and M, respectively. Introducing the expansion (|2.6p to the Navier-Stokes 
equation and projecting it over a similar set of test basis fields leads to an ordinary 
differential equation for the coefficients ai nm , which is solved using a 4th-order semi- 
implicit integration scheme. 

The circumference of the outer cylinder in units of the gap is given by 2wro/d = 
27r/ (1 — rf). Already for r\ — 0.71 this is larger than 21, so that a computational domain 
that goes around the cylinder is very elongated. For some of the highest Reynolds number 
simulations we took advantage of the fact that also correlations in azimuthal direction 
decay and restricted the flow to periodically continued azimuthal domains of size 2tt /n sym 
with risym up to 9. In the spectral expansion (|2.6|) this is achieved by substituting n = 
fi S y m n' and summing over n' instead of n. Then all functions are 27r/n SJ(m -periodic in 
the azimuthal direction and repeat n sam -times to fill the complete circumference. In 
the turbulent case, the azimuthal velocity correlations decay over a short length and 
therefore only an azimuthally shorter sub-domain is required to capture the essential 
flow dynamics. A reference calculation for the case of the outer cylinder at rest showed 
that the difference in the torque was less than 1%. Therefore, a suitable choice of n sym > 1 
substantially decreases the computational costs and gives access to the highest Reynolds 
numbers. 

2.3. Convergence criteria 

Convergence of the calculations is tested using three criteria. 

The first one exploits the radial independence of the current definition (|2.4[) and com- 
pares the values obtained for the inner and outer cylinders. A calculati on is considere d 



conve r ged w hen the two values agree within their fluctuations, see also iMarcusI (|1984[) : 
iDond (i2007h . 

A second convergence criterion requires the coefficients in the Fourier-Chebyshev ex- 
pansion (|2.6p to be sufficiently small. In order to define a measure of the relative drop of 
the coefficients in different physical directions, we normalise the maximal amplitude of the 
highest mode in each direction by the globally strongest mode a max = max; n m \ai nm \, 
i.e. 

d L = max{|a Lnm |, |a_ 

n.m 

a N = max{|a ;A r m |, \ai_ Nm \} / a mwL , (2.7) 
a M = max \a LnM \/a max . 

1,71 



By this definition the amplitudes a^, cln and clm provide a measure for the dynamical 
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amplitude range covered in the corresponding expansion. Additionally, they enable a 
comparison of the quality of approximation in the three directions. For a properly re- 
solved simulation, the values of cll, a at and &m should be of comparable magnitude and 
sufficiently small. 

A third criterion exploits energy balance: The energy injected by the torque and the 
driving of the cylinders at constant velocity must equal, in the statistically stationary 
state and as a n avera ge in time, the volume energy dissipation rate e. As shown in 
lEckhardt efall (|20076h this gives a relation between i and the rescaled torque Nu^, 

i = Re 2 s Nu bJ with e = d 4 v~ 2 ((diUj + d J u i ) 2 ) vt . (2.8) 

Note that while (|2.8p follows as an exact relation from the full Navier-Stokes equation, 
it is not satisfied in under-resolved numerical simulations and can therefore serve as 
another convergence requirement. It was previously employed by iMarcu s (1984) as a test 



of numerical a ccura cy and its meaning for Rayleigh-Benard convection is discussed in 
Stevens et al\ (|2010l ). 



As a first test, we consider a simulation of turbulent Taylor vortex flow for Rei — 
5000 and stationary outer cylinder. We start with a basic resolution of (L, N, M) = 
(46, 70, 30) that shows a comparable resolution in all spectral directions (cf. figureQJb)) as 
a reference and scale the number of modes up and down for different resolutions. Figure [T] 
shows the results for the three convergence criteria mentioned before. First, one observes 
that the quasi-Nusselt number is overestimated by under-resolved simulations, which 
was also noticed for Nusselt numbers in the simulation of Rayleigh-Benard convection 
( Stevens et all l2010h . The three highest resolutions result in approximately the same 



torque which is taken as an estimate of the asymptotic value and indicated by the dashed 
grey line in figure [U a). Torques calculated at the inner and outer cylinder overlap for a 
resolution of (L,N,M) = (36,56,24), but since the approach to the asymptotic values 
is monotonic from above, this agreement is misleading and a sufficient approximation to 
the asymptotic value is only reached with the next step in resolution. 

For resolutions above (L, N, M) = (46, 70, 30), where the torque values does not change 
anymore, the amplitude of the highest mode in each spectral direction drops to 10 -4 , see 
figure Eh). We therefore require that a range in amplitudes of at least 10 4 is covered by 
the spectral modes. 

The third criterion on the relation between torque and energy dissipation turns out to 
be the most demanding: for the resolution (L, N, M) = (46, 70, 30) it is only satisfied to 
within 1% and the error falls of rather slowly with increasing resolution. 

An analogous analysis of the convergence criteria carried out for the case of counter- 
rotating cylinders {Rei = —Re = 2500) resulted in similar findings. Following these 
convergence characteristics, all subsequent simulations conform with three convergence 
criteria: agreement of torque measurements at the inner and outer cylinder to a relative 
deviation of 5 x 10 -3 , resolution of all spectral directions with relative amplitudes of the 
highest modes ~ 10 -4 and fulfilment of the energy balance (|2.8[) to within 10~ 2 . 



3. Effect of vortex sizes on the torque 

Computation of axis-symmetric Taylor vortex flow ( Riecke fc Paap Il986l ) show that 
the torque depends on the axial size of Taylor vortices. Moreover, not only the actual 
torque but also the scaling with Reynolds number va ries for different sizes o f turbulent 
Taylor vortices as previously measured (see figure 2 of ( Lewis k, Swinnevlll999[) ). Since we 
keep the aspect ratio small so as to be able to reach higher Reynolds numbers, we studied 
the aspect ratio dependence for the radius ratio rj — 0.71, once for the case of Taylor 
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(24,36,16) (30,46,20) (36,56,24) (46,70,30) (62, 94, 40) (78, 126, 50) 

(L, TV, M) 

Figure 1. Convergence analysis for a simulation with Ra = 5000 and stationary outer cylinder 
assessed by (a) the torque computed at three radial locations, (b) the amplitudes of the highest 
spectral modes and (c) the deviation from the energy balance (|2.8[1 . In (a) the error bars for 
Nuu,(r a ) denote the confidence interval due to temporal fluctuations. At the cylinder walls these 
confidence intervals are approximately of symbol size and are left out for this reason. The torques 
computed at different radii coincide for the highest four resolutions in (a). 



vortex flow (TVF) and then for turbulent Taylor vortices (TTVF). The procedure used 
in both cases is the following: Simulations at the same Ret were started for various values 
of the aspect ratio T. Then the states that formed were characterised by the number of 
vortices n v identified from flow visualisation. We continued by increasing or decreasing 
r in small steps, with the flow field of the previous T-simulation as initial condition. The 
step size in terms of the axial wavelength of a Taylor vortex pair A z = 2T/n v did not 
exceed 0.1. 

For the first case of TVF, the Reynolds number Ra = 150 was chosen for two reasons: 
It corresponds to Rei/Re c = 1.85, where Re c denotes the critical Reynolds number for 
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Figure 2. The rescaled torque in dependence on the axial wavelength of vortices A z (a) for the 
case of Taylor vortex flow (Rei — 150) and (b) turbulent Taylor vortex flow (Rei = 5000), both 
with stationary outer cylinder. The torques were acquired from states with various numbers of 
vortices. While the error bars in (b) exemplify statistical uncertainties due to fluctuations, Nriu 
was constant in (a) and no error bars are shown. 



TVF, and is thus considerably beyond the bifurcation point. This implies that a broad 
band of Eckhaus stable wave numbers is available and each wave num ber corresponds 
to a s table Taylor vortex state of the corresponding axial wavelength Riecke fc Paapl 
(1986). On the other hand it is below the secondary bifurcation to wavy vorte x flow, 
which occurs for Rei/ Re c > 2.0 according to a numerical analysis for r) = 0.75 ( Jonesl 
19851) . The chosen Rei = 150 is therefore in between the first and second bifurcation 



point and allows for the study of vortices of various sizes. 

As a first step, states with 2, 4 and 6 vortices of the same axial wavelength A 2 = 2 
were simulated in aspect ratios T = 2, 4 and 6, respectively. The first two states were 
continued by increasing T. This resulted in bifurcations going from 2 to 4 and from 4 
to 8 vortices. Note that the latter bifurcation skipped the 6-vortex state and led to an 
8-vortex state at r = 6. These (now four different) states were followed down in T, and 
showed transitions from 8 to 4, from 6 to 4 and from 4 to 2 vortices. The 2-vortex state 
was continued to the lowest realised aspect ratio of T = 1. 

The corresponding rescaled torque values are shown in figure [2ja). Plotted against 
the axial wavelength A z of vortices, the torque values of the different states collapse, as 
expected from the axial periodicity of the DNS code. Nu^ in dependence on the vortex 
size exhibits a broad smooth maximum. The transport of angular velocity is maximised 
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for Az = 1.68, which is smaller th an the critical wavelength A c = 2.0 at the bifurcation 
to TVF (jRecktenwald et aZ.lll993l ). 

The same kind of computations for turbulent Taylor vortex flow at Rei = 5000 (with 
n sy m = 6 to make the simulations feasible) led to similar results, see figure [^b). This 
time, we observed transitions between 2 and 4 vortices and from 6 to 4 vortices. However, 
a transition from 4 to 6 vortices did not occur within the investigated interval up to 
r ^ 6.6. Again, the Nu^ values collapse reasonably well in terms of A z . The torque is 
maximised for X z = 1.93 which means an increase compared to the previous low-i?e; 
value. Furthermore, the functional shape of the torque is markedly different, with a more 
pointed maximum, but smaller relative variations overall. 



4. The angular velocity current 

The global torque is obtained as the spatial average J" over a local quantity, the local 
current j u (r, <p, z,t) (|2.3p . The spatial and temporal fluctuations of the local quantity 
are influenced by turbulent flow-structures and show different characteristics near the 
walls and in the middle of the gap. However, the averages agree as required by the radial 
independence of the mean J". We here analyse local properties of j u for a turbulent 
Taylor vortex flow at Res = 19737 and stationary outer cylinder. 

The spatial variations of the local current j u adjacent to the cylinder walls and on a 
surface at midgap are shown in figure Eta)-(c) for a flow field of turbulent Taylor vortex 
flow. Note that the instantaneous torque results from an average of j u over such cylin- 
drical surfaces. Comparing the section at midgap (figure EJb)) with the arrows in figure 
[5] reveals that the horizontal lines of low and high current j" correspond to the inflow 
and outflow regions, respectively. The spatial j u fluctuations at midgap are dominated 
by these line structures and result in extreme events that are two orders of magnitude 
stronger than at the cylinder walls and than the mean J" = 2.62 x 10 6 . Such extreme 
fluctuation amplitudes have also been seen experimentally, where it has been noted that 



the l ocal current can exceed the mean by a factor of more than 300 (Hu isman et al 
20121) . 



At the origin of the outflow boundary a herringbone-like pattern of low current regions 
can be seen (figure [3f a) ) - This patte rn corresponds to the high- velocity streaks previously 
reported for rj = 0.5 in bond (l2007h : High azimuthal velocity at some distance from the 
wall leads to a small radial derivative of u) at the inner cylinder surface. Between these 
herringbone stripes one can find blobs of extremely high j w current originating from a 
local steepening of the radial angular velocity gradient. 

Similarly, at the outer cylinder (figure EJc)), one recovers V-shaped stripes at the 
origin of the inflow boundary with high j u blobs in between. Note that the V is split 
into one part near the top and one near the bottom boundary in z, but that they are 
connected because of the periodic boundary condition. Here these stripes coincide with 
low-speed streaks and both stripes and blobs are axially somewhat broader than at the 
inner cylinder. In addition, a band of high j u current can be observed at the outer wall 
where the strong axially localised outflow ends and causes a steepening of the w-profilc. 
Interestingly, the counterpart at the inner cylinder wall seems to be absent, which may 
be due to a weaker inflow. 

Averaging the local transport of angular velocity over the corresponding cylindrical 
cross-sections results in the instantaneous torque which still exhibits fluctuations as 
shown in the time-series figure |3jd)-(f). Again the fluctuations in the middle of the 
gap are largest. The mean torque, included as a grey line in the time-series, is obtained 
by additionally performing the time average for the torque signal at each radial location 
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Figure 3. Spatial and temporal fluctuations of the current j u resulting in the torque for 
Res = 19737 and a — 0. Cylindrical cross sections of j u at (a) the inner cylinder, (b) at 
midgap and (c) at the outer cylinder are shown and simulated for one sixth of the annulus. The 
abscissa always covers an angular range of (0, 2ir/6) but is scaled so that the lengths reflects that 
the increase of the cylinder circumference with radius. The corresponding temporal fluctuations 
for the averaged current are shown as Nusselt numbers in (d), (e) and (f). Note that the scale 
in (b) is about a factor 100 larger than in (a) or (c). The mean current is J" = 2.62 x 10 6 , the 
mean Nusselt number Nu u = 15.7 marked by a horizontal line in (d), (e) and (f). 
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Reynolds stress-like current j" s — r 3 v r lj v 

convective current of the laminar profile j" = r 3 v r fit, 

viscous derivative (molecular current) j" d = — v (r 3 d r 0Jv — 2Ef) 
total ^-current = ff s + j" + j" d 

Table f . Local currents (scalar fields) contributing to . The names reflect their physical 

meanings. 



independently. We observe the same mean torque value despite of the different fluctuation 
behaviours. 

4.1. Contributions to the current 

The two contributions of convective and molecular transport in the angular velocity 
current (|2.4[) can be further decomposed when separating the velocity field u into the 
laminar base flow Ub and deviating field v, 

u = ut + v = v r e r + r (fij + w„) e v + v z e z with J7f, = A + B/r 2 . (4-1) 

With this decomposition (|2.3I) becomes 

f = r 3 (v r u) v )+r 3 (v r Q, b ) -vr 3 d r u) v + 2jaB ; (4.2) 

(1) (2) (3) (4) 

now with four contributions: (1) the turbulent convective transport results in a Reynolds 
stress-like contribution; (2) the radial transport of the laminar profile may be non-zero 
locally, but vanishes upon averaging due to the incompressibility condition which implies 
( v r)A(r) = 0; (3) the viscous derivative of the deviating angular velocity; and (4) the 
laminar term which contributes a constant. This motivates the definition of three fields, 
where the last two terms in (|4.2|) are combined, see table [T] As mentioned, the second 
term in (|4.2p contributes to the local fluctuations of the w-current, but has no impact on 
the instantaneous net transport through a surface A(r) and, thus, does not contribute to 
the torque. But it does have a considerable impact on the local fluctuations, as we will 
see in the flow visualisations in the following. 

We analyse the contributions of the different terms for the example of turbulent Taylor 
vortex flow with Reg = 19737 and stationary outer cylinder which was already shown 
in figure |3] For this case the corresponding instantaneous fields are visualised in a r-z- 
plane for a fixed angle ip in figure |4] and EJ First the downstream velocity component 
lo = u v /r is shown with in-plane components represented by arrows (figure 0]), followed 
by a visualisation of the current j u (figure Ufa)). Then the different contributions to 
the current are presented (figure E(b)-(d)). Using this flow case, some properties of the 
distribution of the contributions to j u can be discussed. 

Even for the Reynolds number Res = 19737, an inflow and outflow region is noticeable 
in the instantaneous flow visualisation (arrows in figure [31(a)) indicating the presence of 
turbulent Taylor vortices. The flow shows a transition to the angular velocity prescribed 
by the cylinders in an extremely narrow region adjacent to the inner wall and a broader 
region near the outer wall, see figured Mushroom-like structures of angular velocity are 
pulled out from the wall, especially at the inner cylinder where they are smaller (cf. the 
magnifications in (b) and (c)). In a cylindrical plot, the mushroom s would correspond to 
azimuthally aligned h igh-o; streaks also repo rted bv iDongl (2008 a). Within the scope of 



the w-current theory (jEckhardt et a/Jl20076l ) , which is in analogy to thermal convection, 
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Figure 4. Cross section of instantaneous angular velocity at a fixed angle for Res = 19737 and 
a = with y = (r — ri)/d. (b) and (c) are details of the marked regions in (a) at the outer and 
inner cylinder wall, respectively, with arrows indicating the in-plane velocity components. 



these w-detachments from the boundaries correspond to thermal plumes in Rayleigh- 
Benard flow. 

The flattening of the angular velocity to a medium level in the bulk leads to a negative 
(positive) deviation from its laminar value in the inner (outer) region. Therefore, radial 
outflow implies a negative contribution to the net current near the inner cylinder 
and a positive contribution near the outer cylinder. For the radial inflow, the signs are 
reversed, see figure [5jb). In contrast, the consideration of the convective transport of ilb 
(not contributing to the mean J w ) leads to a current distribution that coincides with 
the radial flow, see figure EJc) . Since ff s generally is much smaller than the convective 
transport of the laminar profile j£, it only weakly influences the complete convective 
current j" a + j and the total current j u , cf. figure [5fa). Compared to j", the viscous 
derivative contribution to the current generally is about two orders of magnitude smaller 
and confined to the cylinder walls to such an extent that areas of maximal j^ d are hard 
to identify in figure [5Jd). This indicates that strong radial w-gradients only occur in ex- 
tremely narrow regions near the walls (cf. figure|3]). Most of the recognisable structures of 
3vd m ^ ne bulk are of negative sign, much weaker in amplitude than the ones expected at 
the walls, and negligible in relation to the convective current contributions. In summary, 
the spatial structures of the total current j w are dominated by the convective contri- 
bution and thus by the occurrence of radial flow (weighted with r 3 f2f>), see figure EJa). 
Consequently, the main part of the local current fluctuations (also in figure |3jb)) origi- 
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Figure 5. Cross sections of local currents at a fixed angle for the flow situation depicted in 
figure 31 Figure (a) shows the angular velocity current, (b), (c) and (d) show the terms that 
contribute to (a) according to table [T] Arrows indicate the in-plane velocity components. Note 
that the scales in (b) and (d) are about 10 times and 100 times, respectively, smaller than the 
one in (c). The local fluctuations exceed the mean J" = 2.62 x 10 6 by two orders of magnitude. 
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nates from the transport term, which, however, does not contribute to the torque at all. 
These large contribut ions may also explain the extreme local observed experimentally 
( Huisman et al. 2012i h Fluctuations of the convective current contributing to the torque 



are much smaller, and only ~ 30 times larger than the mean J w (figure [H^b)). 

4.2. Spatial fluctuations 

The local values of give the local shear stress and their probability density function 
(PDF) near the wall can be compa red to the PDF of the l ocal wall shear stress that was 
me asured at the inne r cylin der bv lLathrop et al. (1992a) and later at both cylinder walls 



bv lLewis fc Swinnev ( 19991 ). Both studies reported log-normal distribution functions in 



these cases. 

In order to compare to th e experimental situation, the data for Res — 81555 and a sta- 



in uiuvi uv vuiii|juii v_< u \j uiiv wi|yvjiim\jiituii uiuumuiuii, uin- uuuu awa j. ^ i_> a t_/ \s t_/ ujaava uj u uuj 

tionary outer cylinder from lLewis fc Swinnevl ( 1999 . fig. 16) were digitised and included 



in Figure [5J where they are shown together with the corresponding PDFs (solid lines) 
computed for the highest Reynolds number achieved in the simulations (Res = 29605) 
and for the outer cylinder at rest. Here, the number of evaluation points was 3/2 times 
the number of modes in azimuthal and axial direction as needed for dealiasing dBovdl 
l200Ct p. 94). One notes that experiment (open circles) and numerical data (continuous 
lines) do not agree, and that the numerical data cover a wider range, including negative 
values. One also notes that while the fluctuations at inner and outer cylinder seem to be 
rather similar for the numerical simulations, they are different for the experimental data, 
with the ones for the outer cylinder being wider. 

To determine whether the differences originate from a .Reg-dcpcndcncc, we show in 
figure [Jj the width (standard deviation a) and the skewness S of the j u fluctuations in 
relation to the mean J u for simulations at various Res- We find that the width of the 
distribution at the inner cylinder varies only slightly with Res- This agrees with the 
experimental observation that the standard deviatio n of the wall shear st ress fluctua- 



tions at the inner cylinder is proportional to its mean ([Lathrop et a/.lll992al ). The width 
of j u fluctuations at the outer cylinder decreases slowly with Res and approaches the 
inner cylinder level (figure [7]) • For completeness also the skewness S of the fluctuation 
distributions is given: While the outer cylinder skewness follows the trend of the corre- 
sponding standard deviation, the one for the inner cylinder starts with a negative value 
and increases with Res- In conclusion, the pronounced width difference between experi- 
mental and computed j" distributions (figure [6]) is unlikely to come from the difference 
in Reynolds numbers. 

Another difference between experiment and numerics is that in the numerical sim- 
ulations we can obtain pointwise values (within the numerical resolution) whereas in 
experiment there are spatial and temporal resolution issues to consider. In order to ex- 
plore the impact of a sp atial averaging, we no te that based on the dimensions of the hot 



film probe as given by ([Lathrop et q/.lll992al ) the experimental data are obtained with 



a spatial resolution (5 V ,5 Z ) in units of the gap width of (4.9 x 10 3 , 18 x 10 3 ). In the 



numerical simulations the grid point spacing gives a resolution of (14 x 10 -3 , 4.6 x 10~ 3 ). 
These resolutions must be evaluated relative to the observation that the axial correla- 
tion length of j u fluctuations is shorter than the azimuthal one, as is evident from the 
smaller axial characteristic size of turbulent structures in the j u visualisations in Figure 
02 Thus, a relatively small 5 Z is especially important in order to resolve turbulent struc- 
tures, and the experimentally achieved resolution would already have to be improved 
at the Reynolds number of the simulations. Moreover, with an inner cylinder rotation 
of about 10 H z and p rovided a frequency response of the hot film probe of 10 kHz, cf. 
Lathrop et all ( 1992a ). the probe moves 17 x 10~ 3 d during one data sample and there- 
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Figure 6. Distribution of the spatial angular velocity current fluctuations obtained at the 
inner (a) and the outer cylinder (b). The circles correspond to measurements for Res = 81555 
l|Lewis fc Swinnevll 19991 ) and the solid lines give numerical computations for Res = 29605. The 
dashed lines indicate the same numerical data additionally averaged over (n v x n z ) — (7 x 14) 
grid points. 
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Figure 7. Standard deviation a x (solid lines) and skewness S x (dashed lines) of j u fluctuations 
(rendered dimensionless with J") at the inner cylinder [x = i) and the outer cylinder [x = o) 
for different Rei and the outer cylinder at rest; here Res = 2/(1 + rf)Rei. 
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fore more than its width 8 V = 4.9 x 10~ 3 , which may lead to an additional azimuthal 
averaging. Finally, because of the three times larger Res, extreme events occurring in 
the experiment are more localised than the ones in the simulation and require a higher 
resolution for detection. These estimates suggest that the axial and possibly azimuthal 
averaging because of the relatively low resolution of the hot film probe leads to a can- 
cellation of extreme j" events in the experimental distribution, and hence to a narrower 
distribution. Moreover, because of the smaller length scales along the inner cylinder the 
effect should be more noticeable there, consistent with the asymmetry between inner and 
outer cylinder in the experimental data. 

Since it is easy to introduce a spatial averaging in the numerical data, we computed 
j u fluctuations averaged over 7 x 14 grid points in the azimuthal x axial direction. 
The PDF's obtained this way are shown as dashed lines in Figure [6] and reproduce 
the experimental data rather well, with almost perfect agreement for the inner cylinder 
and small deviations at the tails and the maximum for the outer cylinder. The poorer 
agreement for the latter distribution may originate from the stronger dependence of 
the standard deviation of j u / ' J u on Res for the outer cylinder (cf. figure [7]). Overall, 
averaging over plausible domains can account for a large part of the difference between 
experiment and simulation. 

4.3. Temporal fluctuations 

After analysing the spatial j w fluctuations we now turn to the temporal aspects of the 
current, after averaging over cylindrical surfaces A(r). We focus on three locations, the 
inner cylinder, the outer cylinder and the middle of the gap, where the convective contri- 
bution in (|2.3[) dominates. We keep the shear Reynolds number Res — 19737 fixed and 
vary the ratio of the rotation rates, taking a — —0.4, 0.0, 0.4, and 2.0 as representative 
examples. The PDF's for the fluctuations around the mean, normalised by the mean, are 
shown in Figure |8] 

As for the spatial j u distribution (figure |3|) , the fluctuations are again significantly 
larger at midgap than near the inner or outer cylinder. For the first three values of a the 
fluctuations at inner and outer cylinder agree and the fluctuations are symmetric. For 
a = 2 the distributions at inner and outer cylinder arc different and the ones at midgap 
and outer cylinder are strongly asymmetric. We also note that the fluctuations at the 
inner cylinder show little variation with global variation. The fluctuations at midgap are 
strongest for outer cylinder at rest and are smaller for moderate co- and counter-rotation. 
So what happens at a = 2? 

A first indication of a change in flow dynamics is provided by the radial profiles of 
angular velocity for the four rotation states discussed before, see figure [9] Narrow gradi- 
ent regions adjacent to both walls are clearly discernible for co- and moderate counter- 
rotation and stationary outer cylinder. In between the angular velocity exhibits a slight 
decline. The small regions near the walls indicate thin boundary layers. In contrast, the 
angular velocity for strong counter-rotation exhibits a smooth transition from uj to its 
bulk level. One may suspect that this slow transition is due to a low turbulence level 
resulting in a thick boundary layer near the outer wall. This would suggest a connec- 
tion to the known fact that for counter-rotat ion only a region nea r the inner cylinder 
is unstable according to Rayleigh's criterion (jChandrasekhar 1961). The stable region 



adjacent to the outer cylinder increases with increasing counter-rotation. Although this 
stability analysis only holds for the laminar flow, one expects that the turbulence is less 
strongly driven near the outer cylinder. On the other hand, a laminar flow as such would 
be incompatible with the fact that the current (|2.4[) has to be the same, independent of 
the radial position along which it is calculated. The strong fluctuations observed in the 
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Figure 8. Temporal fluctuations of the angular velocity current for Res = 19737 and four 

rotation ratios a computed at three radial locations: , inner cylinder; , outer cylinder; 

, midgap. Except for a = 2.00 the PDFs of the fluctuations at the inner and outer cylinder 

lie almost on top of each other. Note that the scale in (d) is slightly increased. 
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Figure 9. Profiles of the mean angular velocity for the same shear Res = 19737 and four 
different rotation ratios a with y — (r — r»)/d. 
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Table 2. Realised shear Reynolds numbers Res with range of rotation ratio and employed 
resolution. For the outer cylinder at rest or alternatively the rotation ratio specified in the 
second column, the last column gives the quasi-Nusselt number computed from the dissipation 
rate in relation to Nu^ which is suitable as convergence criterion. 



current at the outer cylinder (cf. figure EJd)) are thus needed to compensate in the time 
average the lower than average torque contribution in the less turbulent phase. These two 
flow states a re then consist e nt with the bi-modal behaviour in the outer region recently 
observed bv Ivan Gils et all (12011 6l) for Res ~ 10 6 . Similarly, at much lower Reynolds 
numbers, ICoughlin fc Marcus! (|l99a) described oscillatory changes between laminar flow 
and turbulence in counter-rotating Taylor- Couette flow. 



5. Variations of torque with rotation rates 

Following the local characteristics we now turn to the variation of the global torque 
with rotation rates. More precisely, we study how Nu u , the torque expressed in units of its 
laminar value, depends on the external driving, specifically on the shear Reynolds number 
(Res) and on the global rotation s t ate (a = — uj /u}j). We w ill use the experimental 
observation of ( van Gils et aLll2011al : lPaoletti fc LathroijEml) that the global variation 
can be represented by a scaling in Res multiplied by an amplitude that is a function of 
a only. 

For this purpose we realised 72 simulations at nine shear Reynolds numbers ranging 
from Res — 1155 to Res = 29605 that are spaced by factors of 1.5. The lower Reynolds 
number bound was chosen to ensure reasonably developed secondary flow, whereas the 
upper bound was dictated by computational restrictions. For each shear Reynolds num- 
ber, eight different global rotation states with the focus on the regime of counter-rotation 
were selected. We include here two cases for co-rotation (a = —0.4 and a = —0.2), one 
with a stationary outer cylinder (a = 0), two for moderate counter-rotation (a = 0.2 
and a — 0.4), one for equal but opposite rotation Re = —Rei (a = rj = 0.71) and two 
for high counter-rotation (a — 1.0 and a = 2.0). For the highest Res only one sixth or 
one ninth of the domain was simulated in the azimuthal direction (cf. table [2]). Calcu- 
lations for a situation with stationary outer cylinder show that the differences to a full 
simulation with n sym = 1 are below 1%, so that the reduction is justified. The resolution 
(L, N, M) employed for each Reynolds number fulfils the exact computation of the global 
dissipation rate (|2.8|) within a relative error of 0.7% at most, see the last column in table 

El 

For these flow cases the torques were computed at the inner cylinder and presented 
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Figure 10. Torque dependence on the shear and the global rotation: (a) Compensated plot of 
Nuu, vs. rotation ratio a for various Res. Flow states represented by symbols of the same shape 
feature the same shear, (b) Scaling exponents determined for each rotation ratio independently 
considering torques for Res = 1.32 x 10 4 , Re s = 1.97 x 10 4 and Re s = 2.96 x 10 4 . 



by means of Nu u compensated with Re s ' 76 in figure [TUf a) . The exponent . 76 is t aken 
from the effective torque scaling Nu^ ~ Re° s ' 76 found by Ivan Gils et al. (2011a) for 
much higher Reynolds numbers between 3 x 10 5 and 2 x 10 6 . Each data point on a 
line in figure HOT a) corresponds to a different global rotation state exhibiting the same 
shear. The changes with increasing Res indicate that no global power law Nuqj ~ Re Z 
can be identified, as noted previously by ( Lathrop et al. 1992a ; Lewis fc Swinnev 19991 ) ■ 
Nevertheless, the numerically determined effective torque scaling ex ponent, 7 = 0.67, 
for a s tationary outer cylinder agrees well with 7 = 0.69 measured bv lLewis fc Swinnev 
(1999) at Res ~ 1.5 x 10 4 . In addition to the variation with Res, the local scaling 
exponent 7 depends on the rotation ratio, see figure [TU^b). The exponent clearly decreases 
from moderate counter-rotation towards co-rotation, i.e. with decreasing rotation ratio. 
A similar tre nd is also discernible in effective exponents from high Reynolds number 
experiments (jvan Gils et a/Jl2011M ). However, the authors note that both a trend in the 
exponents and a statistical scatter around an a-independent exponent are within the 
resolution of their experimental data. Presumably, the variations of 7 in figure HUT b) 
are larger than in the experiment because of the relatively low Res ~ 10 4 , where the 
flow is not yet turbulence-dominated. Note, however, that the maximal scaling exponent, 
7 = 0.77, in our simul ations for o = . 4 alrea dy agrees well with the ultimate exponent, 
7 = 0.76, measured by van Gils et al. ( 2011a| ). 

Because of the complicated dependence of the torque on the shear and global rotation, 
a crossover in the rotation ratio that maximises Nu^ is observed: Initially, the transport 
contributing to the torque is most effective in a situation with stationary outer cylin- 
der. However as the Reynolds number increases, the torque maximum becomes flatter 
and shifts towards moderate counter-rotation for Res ^ 2 x 10 4 . In terms of the shear 
Reynolds number, this shift of the maximum towards a ~ 0.4 can be seen as a pre- 



cursor for the torque maximisation at a — 0.368 for 77 = 0.716 ( van Gils et al. 201161 ) 
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Figure 11. Comparison of torque betw e en the numerical simulation (rj = 0.7 1, Res = 2.96xl0 4 ) 
and th e experiments bv lvan Gils et ali (|2011al ) (Twente, r\ = 0.716) and bv lPaoletti fc Lathropl 
(|201lT ) (Maryland, r\ = 0.724). (a) shows rescaled torques. Note that in the Twente experiment 
a rotation ratio-independent effective scaling Nu^ ~ 7?eg 76 was found and the a-dependent 
prefactor of this scaling law is plotted here, (b) shows relative torques for both experiments and 
the simulation. The legend also applies to (a). 



and at a = 0.333 for r\ = 0.7245 ( Paoletti fc Lathropl 2011 ) as observed for even higher 
Res- However, both experimental studies did not report a shift of the torque maximum, 
presumably because the Reynolds number range they explore was much higher. 

Observing the maximisation of the transport for approximately the same global rota- 
tion state, it is tempting to directly compare the functional dependence on the rotation 
ratio between the simulation at Res — 3 x 10 4 and the experiment even though they 
are separated by one to two orders of magnitude in Reynolds numbers. Appealing to the 
universal scaling found experimentally and the small change between the two highest 
Res in figure fTCT a). it seems appropriate to compare the rescaled Nu^Reg ' 76 . For this 
type of comparison we assume that the mean dependence on Res is divided out by the 
rescaling. 

In the experimentally found universal torque scaling ( van Gils et al. 2011al ). only a 
rotation ratio-dependent proportionality factor in the power law absorbs the complete 
influence of the global rotation on the transport. This factor is given together with the 
corresponding calculated iVu w for the highest Res in figure [TTT a). Qualitatively, the 
functional dependence on the rotation ratio agrees between simulation and experiment 
and the tails compare well even quantitatively. However, the torque maximum is not as 
pronounced in the simulation as in the experiment. This may be due to the relatively low 
Reynolds numbers accessible by simulations which is supported by the variation of the 
local scaling exponent in figure [TOTb) . Its maximum for a = 0.4 suggests that the torque 
maximum will become more pronounced in simulations for even higher Res- 

A way to compare torque without the assu mption of some scaling use s the torque of 
the inner cylinder at rest as a reference value. IPaoletti fc Lathropl (|201l[ ) observed that 
this torque ratio is almost independent of the shear Reynolds number and fitted the 
dependence on the global rotation by linear functions of Ro^ 1 = u! /(uii — lo ) within 
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four regions. Their functions for region 3 and 4 ( Paoletti k, Lathropl 2011 , eq. (5), (6)) 
are shown in figure [TlT b) together with data from the Twente experiment and from 
simulations. Except for the height of the iVu^-maximum the computed torque values also 
agree with these empirical functions. Additionally, the observations of both experiments 
agree well with each other with a small deviation for strong counter-rotation. Note that 
this comparison only analyses the functional dependence on the global rotation based on 
relative torques. Therefore, possible differences in the absolute Nu^ values can not be 
identified in figure [TlTb). 



6. Relation between torque and boundary- layer thicknesses 

Boundary layers play a key role in the analysis and description of global transport 
properties, and the numerical simulations presented here provide a means to compare 
different definitions of boundary- layer thickness, and to pick the one that is most suitable 
for the development of the theory. The observation that the torque shows no universal 
scaling with Res in the investigated regime and, furthermore, depends on the mean 
rotation indicates that also the boundary-layer thickness will have a more complicated 
dependence than an effective power law in Re s with a rotation-dependent prefactor, as 



obser ved experimentally for much higher Res (|van Gils et aLll2011atlPaoletti fc Lathrop 



20111) . 

Thin layers near the cylinders, where the angular velocity changes from the wall to the 
bulk level, develop as a result of the applied shearing and are already visible in figure 
[9] These angular velocity profiles illustrate boundary layers (BL) due to azimuthal flow 
over the cylinder surface in the rotating frame with the corresponding cylinder at rest. 
Similarly, the axial flo w over the cylinder surf ace is assumed to cause a second BL of an 
independent thickness (jEckhardt et Since the angular velocity profile directly 



enters the torque equation (|2.4p . only the first BL type will be analysed here. 

Common measures of the BL thickness for the flow over a flat plate include the 99% 
velocity thi ckness and displa cement thickness and require a free stream velocity Woo to be 



computed (jSchlichtina 119821 ). However, such an unambiguous free stream velocity does 
not exist in the TC flow because angular velocity profiles exhibit different slopes in the 
bulk depending on global rotation, cf. figure [HI Another measure consists in the definitio n 



of a slope thickness based on the angular velocity, as proposed bv lEckhardt et al. ( 20076 ) 



For this purpose one relates the slope of the w-profile at the boundaries to a drop in 
angular velocity from the surface to a representative bulk value, \lo x — uj\. This leaves the 

(x) 

boundary-layer thickness X sl ' as an unknown, 



d r (u) A . t k = {x ~' 5 x = i,o. (6.1) 

Kl 



Here Q denotes a bulk or centre level of the angular velocity as a point of reference. The 
result for the case of co-rotation (a = —0.40) with the definitions u> = (uj(r a )) At and 
J"o = ( r i +r )/2 shows that this overestimates the size of the region one would intuitively 
identify as BLs, see the dashed lines in figure Q21 The obvious reason lies in the fact that 
the reference value ui is estimated far from the boundary region, thereby giving a larger 
drop in angular velocity, and this has to be compensated by a correspondingly larger 
value in A s ;. 

To overcome this problem, the definition of an "effective boundary-layer thickness" 



Torque in Taylor- Couette Flow 



21 




0.6 ' — "- 1 1 ' 1 1 '—— 1 

0.0 0.2 0.4 0.6 0.8 1.0 

V 

Figure 12. Mean angular velocity profile for co-rotation (a = —0.40) and Res = 19737 to exem- 
plify the definition of boundary-layer thicknesses. The dashed lines indicate the slope thickness 
A s ;, while the dotted lines indicate the effective boundary-layer thickness X v d for the inner and 
outer cylinder BL, respectively. 



with x — i,o and radial locations, r^\ corresponding to 



^(rjW) = min{J^} + 0.1 (> - min{J- } 

r \ r 



(6.2) 



can be employed (|Donell20086h . The subscript vd indicates that it is based on the tran- 
sition between viscous derivative Jy d (r) = (j^AM t anc ^ convective contribution to the 
current. It identifies the distance from the wall where the viscous derivative term exceeds 
its minimal value by 10%. Near this point the slope of the w-profile approaches its bulk 
value, but it does not necessarily vanish (cf. figure [12|) . The definition (16.21) is intrinsi- 



(Hajsjr 
/ lEckhard 



cally cl ose to the flow region description in the torque scaling theory by lEckhardt et 
(|2007a ). The measure X v d corresponds to the BL regions one would expect from the 



shape of the w-profile for a = —0.40, see dotted lines in figure IT21 Note that the relation 
And < for a — —0.40 is inverted for other rotation ratios. 

So far, we introduced two different measures for angular velocity BL thicknesses which 
characterise the turbulent flow situation. Like the torque, these thicknesses vary with the 
externally applied shear and mean rotation. In the following, we analyse to which extent 
the variation of both BL thickness definitions agrees with the functional dependence of 
the torque. More precisely, we ask whether the knowledge of BL thicknesses suffices to 
predict the torque. 

In order to study this issue, we employ the relation 



= a 



A(«) 



(6.3) 



with a weighted average w-BL thickness and a = ((1 + if) / (2 y 7 ?;)) 4 , which results 
fro m adding the two ex pressions (|2.4[) at the boundaries and using equation (|6.ip . 



sec 



Eckhardt et ali (|2007M ). Furthermore, one has to assume that the inner and outer 



boundary-layer thickness are related by 
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Figure 13. Test of relation (|6.3p between torque and average BL thickness A' a ' calculated from 
the slope BL thickness A 3 ; in (a) and from the effective BL thickness X v a in (b). The grey line 
indicates the value a expected from theory. For symbols in (b) see legend in (a). 

which defines the average thickness A^ a ^. Relation (|6.4j) originates from the observation 
that the constant current J u together with (|6.1|) results in 

„3 |, , -I „3 

(6.5) 



A(°) |w -w| AW 



In consequence, the assumption ()6.4|) consists in the approximation that the ratio of 
angular velocity drops \u>i — Q\/\lj — u>\ is close to unity. Here the weighted average BL 
thickness is defined for both BL thickness measures as 



(6.6) 



consistent with relation (16.41) . According to (|6.3[) . the product 2Nu u) X^ / d should be 
constant and close in value to a. 

We first compute X^ using the slope thickness X s i with a properly chosen bulk level 
of the angular velocity u> — {uo{r a )) At . We find that the deviation from expected relation 
(|6.3p is not larger than 11%, see figure [TBT a). This high degree of conformance can 
be expected since the BL thickness measure X s i was defined to be consistent with the 
description of the current J" and requires the knowledge of d r (u) A t at the boundaries. 
The latter, however, can be directly used to calculate Nu u . Nonetheless, the finding is 
not trivial because we needed the assumed relation (|6.4[) which postulates that the ratio 
of angular velocity drops to the bulk level w equals unity. 

The situation drastically changes when Nu^ is instead related to a weighted average 
of the effective BL thickness X v d- The product 2Nu L0 X^ jd shows a transient phase with 
increasing Res, especially for the cases of strong counter-rotation, see figure H3TbV Then, 
the values vary less with Res, but deviate from a and differ for the various rotation ratios. 
A maximal proportionality coefficient is reached for moderate counter-rotation. The large 
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deviations from a suggest that the BL thickness measure X v d is not consistent with the 
characteristics of the current resulting in the torque. 

In summary, the torque can be derived from BL thicknesses in the full parameter 
range considered here. However, the connection only works for a particular measure 
of the BL thickness, namely the slope thickness X s i that is combined in a weighted 
average (|6.6p . Therefore, the knowledge of an angular velocity profile suffices to predict 
the torque. Although the definition of \ v d is based on the current J u , it fails to explain 
the dependence of the torque on global rotation. 



7. Conclusions 

Direct numerical simulations with accurate torque calculations were performed for 
Taylor- Couette flow and achieved Reynolds numbers Res ~ 3 x 10 4 not previously re- 
ported in literature. These simulations enabled comparisons to experimental investiga- 
tions of turbulent Taylor-Couette flow at much higher Reynolds numbers. 

These direct numerical simulations show that Taylor-Couette turbulence in domains 
that carry only one pair of rolls and that cover only part of the circumference can give 
insights into the transition from the vortex dominated flows at lower Reynolds numbers 
(up to about Res ~ 10 ) and the fully developed turbulence, which begins to show up 
near Res = 3x 10 4 , the highest shear Reynolds numbers realised here. While the aspect 
ratio is rather important in the Taylor- vortex regime, it seems to become less significant 
in the turbulent cases. 

The statistical analysis of the local current j w showed that the fluctuations are wider 
than experiments suggested and they reach to negative values. The reduction in width of 
the fluctuations could be explained by an averaging over space: this suggests that future 
measures of these fluctuations require smaller probes. The fact that the fluctuations can 
become negative indicates that locally in space and time there can be fluctuations that do 
not take energy out of the mean flow but actually deposit energy there. This is consistent 
with observ ations on plane shear flow whe re a connection to fluid structures could be 
established ( Schumacher fc Eckharddlibol . 

The observation of the changes in mean profile with strong counter-rotati on extend 
the observations of fluctuating boundary layers from lCoughlin fc Marcus! (|l996l ) to higher 
Reynolds number. The onset of this phenomenon and the thickness of the affected region 
require additional study that will be presented elsewhere. 

We have also identified a definition of boundary-layer thickness that is compatible 
with the global torque rela tion and that can be use d in the corresponding analysis of the 
global transport relations Eckhardt et all ( 20076( ). The definition is more complicated 
than the usual definitions because the reference value in the bulk is not a constant, 
independent of the radius. The variation of the mean profile in the middle with radius 
is ano ther remarkable observation, which, however, could also be gleaned from IWendtl 
(1933). The origin of this variation, which is a difference between Rayleigh-Benard and 
Taylor-Couette flow, remains unclear. 

The comparison to experiments was limited by a large gap in the Reynolds numbers: 
The simulations are limited to Res near several 10 , whereas experiments begin at several 
10 5 . Nevertheless, several features like the development of a torque maximum for counter- 
rotation could be studied. Experimental investigations of the flow at moderate Reynolds 
numbers would be of great value to draw direct comparisons between experiments and 
simulations in the turbulent case. 

The torque maximum, the boundary-layer variations and the fluctuation statistics 
investigated here suggest that studies of turbulent Taylor-Coutte flow may be just as 
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rewarding as the analysis of the many bifurcations that lead from the laminar Couctte 
flow to the turbulent flow state. 
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